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ABSTRACT 

In this letter, we present the resuhs of our study of galaxy-galaxy lensing in massive cluster-lenses 
spanning z = 0.17 to 0.58, utilizing high-quality archival Hubble Space Telescope (HST) data. Local 
anisotropics in the shear maps are assumed to arise from dark matter substructure within these clusters. 
Associating the substructure with bright early-type cluster galaxies, we quantify the properties of typical 
L* cluster members in a statistical fashion. The fraction of total mass associated with individual galaxies 
within the inner regions of these clusters ranges from 10-20% implying that the bulk of the dark matter 
in massive lensing clusters is smoothly distributed. Looking at the properties of the cluster galaxies, 
we find strong evidence (> 3-(t significance) that a fiducial early- type L* galaxy in these clusters has 
a mass distribution that is tidally truncated compared to equivalent luminosity galaxies in the field. 
In fact, we exclude field galaxy scale dark halos for these cluster early- types at > 10-cr significance. 
We compare the tidal radii obtained from this lensing analysis with the central density of the cluster 
potentials and find a correlation which is in excellent agreement with theoretical expectations of tidal 
truncation: log[rt*] oc (— 0.6 ± 0.2) log[po]- 

Subject headings: gravitational lensing, galaxies: fundamental parameters, halos, methods: numerical 



1. INTRODUCTION 

Galaxy-galaxy lensing provides a powerful tool to sta- 
tistically measure the mass and the details of the mass dis- 
tribution for field galaxies (Tyson et al. 1984; Brainerd et 
al. 1996). These studies confirm the existence of massive 
dark matter halos around typical field galaxies, extend- 
ing to beyond 100 kpc^ (Brainerd et al. 1996; Fischer et al. 
2000; Smith et al. 2001b; McKay et al. 2001). The same 
technique can be modified and implemented within clus- 
ters to constrain the masses of cluster galaxies (Natarajan 
& Kneib 1997, NK97; Geiger & Schneider 1998). Suc- 
cessful application of the same to the rich, lensing cluster 
AC 114 at z = 0.31, suggests that the average M/L ra- 
tio and spatial extents of the dark matter halos associated 
with early-type galaxies in such dense environments may 
differ from those of comparable luminosity field galaxies 
(Natarajan et al. 1998, NKSE98). 

The technique applied by NKSE98 quantifies the local 
weak distortions in the observed shear field of massive 
cluster-lenses, as perturbations arising from the massive 
halos of cluster galaxies (for details see NK97). By as- 
sociating these perturbations with bright early-type clus- 
ter members, the relative mass fraction in their halos 
is constrained using a combined x^-maximum likelihood 
method. The strength of this approach is the simultane- 
ous use of constraints from the observed strong and weak 
lensing features. 

The fractional mass in clusters associated with indi- 
vidual galaxy halos has important consequences for the 



frequency and nature of interactions (Moore et al. 1996; 
Ghigna et al. 1998; Okamato & Habe 1999; Merritt 1983). 
The theoretical expectation is that the global tidal field 
of a massive, dense cluster potential well should be strong 
enough to truncate the dark matter halos of galaxies that 
traverse the cluster core. In this letter, we test this expec- 
tation using well calibrated mass models for rich clusters 
at z ^ 0.17-0.58 that utilize the observed strong lensing 
features - positions, magnitudes, geometry of multiple im- 
ages and measured spectroscopic redshifts as well as the 
shear field. 

2. GALAXY-GALAXY LENSING IN CLUSTERS 

The clusters in our study are modeled as a super- 
position of a smooth large-scale potential and smaller scale 
potentials that are associated with bright early-type clus- 
ter members: 0tot — 4>c\us+'^i 4>Pii where 0cius is the poten- 
tial of the smooth component and are the potentials of 
the perturbers (galaxy halos) . The resultant reduced shear 
g, ^, g, and can in turn be decomposed into contributions 
from the smooth clump and perturbers. 

To quantify the lensing distortion induced by ^tot, the 
smooth cluster piece and individual galaxy-scale halos are 
modeled sclf-similarly using the pseudo-isothermal ellipti- 
cal (PIEMD) surface density profile, S(i?) derived by Kas- 
siola & Kovner (1993), with a core-radius ro, a truncation 
radius rt S> and an ellipticity e. These parameters 
are tuned for both the smooth component and the per- 
turbers to obtain mass distributions on the relevant scales 



^We adopt h=Ho/100kms ^ Mpc ^=0.5 and Qo = 0.5, f2o = 1- Our results however, are not sensitive to values of the cosmological 
parameters. 

^The reduced shear g, a complex number, is defined in terms of the convergence ft and the shear 7, as g = 7/(1 — k), can be directly related 
to the measured shape parameter t = 2(;/(l — g*g) 
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(for details sec §2.2 of NK97). Note that rt characterizes 
the scale over which the local potential dominates. Addi- 
tionally, assuming light traces mass, the light distribution 
of the early-type cluster members is related to the mass 
model via a set of the usual adopted scaling laws motivated 
by observations (see Brainerd et al. 1996): 

(To = a*{^)i; ro = ro*{^)^; n = rt*{^)^ . (1) 

.L/^ 

The effect of these assumed scaling laws on the maximum- 
likelihood analysis has been explored in detail using simu- 
lations in NK97. The total mass M and the total mass-to- 
light ratio M/ L of cluster galaxies then scales as follows: 

M oc C7»rt*(— ); — oc a^n*. (2) 

The ellipticity of the mass clumps is set to the observed 
ellipticity of each cluster galaxy. Since, we are sensitive to 
the induced local shear arising due to the integrated mass 
within ri for the chimps, we are not sensitive to details of 
the assumed radial fall-off and the assumed inner slope of 
the profile in the galaxy model. 

The measured local shear signal in the cluster at the lo- 
cation of the i-th perturber is boosted in value due to two 
effects (i) the contribution from the smooth cluster com- 
ponent - 7c in the numerator, and (ii) the (kc + Kp.) term 
in the denominator, which is positive and non- negligible 
in the central cluster region. By taking into account a 
smooth cluster model it is possible to reduce the error in 
the measured shear from the cluster members: 

where gi,gs and Qc are the shapes of the lensed image (I), 
unlensed source (S) and the shear of the smooth cluster 
piece; ^l^rs) ^ '^^ width of the intrinsic elhpticity 
distribution of the sources, and Ni,g the number of back- 
ground galaxies averaged over. Therefore, knowledge of a 
well calibrated strong lensing model makes the study of 
galaxy-galaxy lensing in clusters very viable. 

A maximum-likelihood method is used to obtain sig- 
nificance bounds on fiducial parameters that characterize 
a typical L* halo in the cluster. We have extended the 
formalism developed in NK97 to include the strong lens- 
ing data for the inner regions of the clusters, these are 
used to obtain the best model fit followed by a likeli- 
hood method that incorporates the constraints from the 
smoothed shear field. 

The likelihood function of the estimated probability dis- 
tribution of the source ellipticities is maximized for a set 
of model parameters, given a functional form of the in- 
trinsic ellipticity distribution measured for faint galaxies. 
The entire inversion procedure to solve the lensing equa- 
tion is performed for each cluster within the LENSTOOL 
utilities'*. Using a well-determined 'strong lensing' model 
for the inner-regions of the clusters along with the aver- 
aged shear field and assuming a known functional form for 
p{ts) from the field, the likelihood for a guessed model is 
computed using: ts^ = Tobsj " Ipi ~ Ic for each faint 

fully analytical ray-tracing routine developed by Kneib (1993) 



galaxy j where " 7^. is the sum of the shear contribu- 
tion from Nc perturbers, and 7c is the shear induced by 
the smooth component. 

£(<7*,n,) = n7-'KTs,). (3) 

We compute C assigning the median redshift correspond- 
ing to the observed source magnitude for each arclet (de- 
tails of this procedure can be found in NK97). The best 
fitting model parameters are then obtained by maximizing 
the log-likelihood function I with respect to the chosen set 
of model parameters (tT*,rtH.) and the cluster parameters 
are simultaneously matched in an iterative way. 

3. THE HST CLUSTER-LENSES 

For our analysis we select clusters at 2; > 0.1 for which 
deep, high-quality HST imaging is available and which 
contain spectroscopically-confirmed multiply-imaged high 

redshift galaxies. These lensed features are essential to 
construct a detailed mass distribution for the cluster cores 
(e.g. Kneib et al. 1996; Smith et al. 2001a), while the exis- 
tence of spectroscopic redshifts allows us to calibrate these 
mass distributions onto an absolute scale. This selection 
yields five clusters with redshifts spanning z — 0.17-0.58 
for our analysis: A 2218, A 2390, CI 2244-02, CI 0024-^16, 
and CI 0054— 27. In addition to these five clusters, we 
also include our previous analysis of the z = 0.31 cluster 
AC 114 (NKSE98). Details of the cluster properties are 
given in Table 1 clearly these clusters do not constitute 
a well-defined sample, for example their X-ray luminosi- 
ties span an order of magnitude and their central mass 
densities show a similarly large dispersion. It is this latter 
property which is of most interest for our analysis - and 
the large range spanned by the sample therefore provides 
a powerful test of the variation in characteristics of galaxy 
halos with local environment. 

The five new clusters were observed with WFPC2, four 
through the F814W (/) passband, with A 2218 imaged in 
F702W (R), the total exposure times are listed in Table 1. 
These typically comprise stacks of multiple single-orbit 
(2.1-2.7ks) exposures, each of which is spatially offset by 
an integer number of WFC pixels to allow the removal of 
cosmic-rays and other artifacts. The data were pipeline 
processed by STScI and combined using standard stsdas 
and IRAF routines. For more details of the reduction see: 
A2218, Kneib et al. (1996); A 2390, Pello et al. (1999); 
CI 0024+16 and CI 0054-27, Small et al. (1997). These 
fields reach typical 5-(j detection limits for point sources 
of / ~ 26 or i? ~ 26.5. 

Object catalogs and attendent shape information for 
faint galaxies are obtained using the S Extractor pack- 
age (Berlin & Arnouts 1996) and detection criteria of 
12 WFC pixels above the /Xfl or fix = 25.0 mag arcsec"^ 
isophote after convolution with a 0.3" diameter top-hat 
filter. Selection of the background galaxies in these im- 
ages employs a simple magnitude cut as previously used 
by NKSE98, either R = 23-26 or / = 22.5-25.5 - this 
provides catalogs of typically ~ 300-400 faint galaxies in 
each field. The magnitude limits employed should select 
field galaxies at median redshifts of z ^ 1, well beyond 
the clusters. Visual morphological classifications of the 
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brighter galaxies in these fields, used to select the early- 
type cluster galaxies, comes from Small et al. (1997, 2001). 

In addition to the strong lens model of AC 114 em- 
ployed by NKSE98, lens models for three of the other 
clusters have been published previously: A 2218, Kneib et 
al. (1996); A 2390, Pello et al. (1999);' CI 0024+16, Smail 
et al. (1996). For the remaining two lenses CI 2244-02 
and CI 0054— 27, we construct well constrained mass mod- 
els. The mass distribution in CI 2244— 02 has signifi- 
cant substructure, with the main clump centered around 
the brightest cluster galaxy with a velocity dispersion of 
600 ± 80 kms-\ an ellipticity of 0.17 ± 0.1 and the sec- 
ondary clump offset about 20 arcseconds away with a ve- 
locity dispersion of 300 ± 40 kms~^ and an ellipticity of 
0.1 ± 0.05. The cluster CI 0054-27 by contrast has a 
smooth, but extended mass distribution (with a cut-off 
radius of about 900±60 kpc) with a central velocity disper- 
sion of 1100 ± 100 kms"^ and an ellipticity of 0.1 ± 0.05. 

4. RESULTS 

Before discussing the results, we describe the various 
tests we have employed to confirm their reliability. In 
order to check the authenticity of the signal the follow- 
ing null tests were performed in the likelihood analysis for 
every cluster: (i) randomizing the positions of the back- 
ground galaxies, (ii) randomizing the orientations of the 
background galaxies, and (iii) scrambling the positions of 
cluster galaxies. No significant maxima were obtained in 
the likelihood prescription for any of these tests in any of 
the clusters. To quantify the errors on our derived pa- 
rameters we note that the principal sources of error in the 
above analysis are (i) shot noise - we are inherently limited 
by the finite number of sources sampled within a few tidal 
radii of each lensing cluster galaxy (in the present analysis 
we typically have 40 cluster members and ~ 400 back- 
ground sources); (ii) the spread in the intrinsic ellipticity 
distribution of the source population; (iii) the unknown 
source redshifts and (iv) observational errors arising from 
uncertainties in the measurement of ellipticities from the 
images for the faintest objects. From simulations where 
the detailed mass distribution of clusters is known, we have 
calibrated the sources of error and find that of the above- 
mentioned four factors, shot noise (accounts typically for 
~ 50% of the incurred error) and the unknown redshifts 
for individual background galaxies (contributes ~ 0-30% 
depending on the redshift of the cluster) dominate the er- 
ror budget. The remaining ^ 20% of the error budget 
arises due to inaccuracies in the measurement of shapes. 

We show in Fig. 1 the likelihood contours for the galaxy 
perturber models of each of the five clusters. In all cases 
we detect an unambiguous galaxy-galaxy lensing signal at 
the >3-(T level - confirming the existence of truncated dark 
halos associated with early-type galaxies in clusters. The 
likelihood analysis yields best-fit model parameters: a* 
the central velocity dispersion and truncation radius 
for a typical L* cluster member, these values are listed in 
Table 1.^ The mass-to-light ratios quoted here take pas- 
sive evolution of the stellar content of elliptical/ SO galaxies 
into account modeled using the stellar population synthe- 
sis models of Bruzual & Chariot (1993). 

*The mass obtained for a typical bright cluster galaxy by Tyson et al. (1998) using only strong lensing constraints inside the Einstein radius 
of the cluster CI 0024+16 is consistent with our results. 



5. DISCUSSION AND CONCLUSIONS 

We have statistically extracted characteristic parame- 
ters for typical L* cluster galaxies that inhabit massive, 
dense lensing cluster-lenses ranging in redshift from 0.17- 
0.58. This has been achieved by combining strong and 
weak lensing HST observations in conjunction with an as- 
sumed parametric mass model. We find that the inferred 
mass distribution of a fiducial L* is extremely compact, 
although the inferred 's lie well outside the optical radii 
and correspond to roughly between 5-10 i?e- Our analy- 
sis also shows that the halos of individual cluster galaxies 
contribute at most 10 20% of the total mass of the clus- 
ter within the central 1 Mpc, covered by the HST WFPC2 
imaging using the results of our likelihood analysis along- 
with the best-fit paramters that characterize the smooth 
clump. Therefore, in the inner regions of these clusters the 
bulk of the dark matter is in fact smoothly distributed. 
Similar lensing studies of field galax;ies, e.g. Wilson et al. 
(2001), typically find a non-zero signal for the radially av- 
eraged stacked tangential shear out to 200 kpc. In contrast 
our study of the halos of galaxies in clusters detects a fi- 
nite rj , which we attribute to the tidal truncation induced 
by the motion of these cluster galaxies inside the poten- 
tial well. From the contours in the likelihood plots, the 
presence of field galaxy scale dark halos can, in fact, be 
excluded at > lO-cr significance. 

The clusters we study here are all rich systems span- 
ning a range in central density, which may explain why 
the best-fit values of obtained vary by a factor of 2-3. 
To test this suggestion we plot in Fig. 2 the variation of 
the central density of the cluster dark matter with r*^ / a* 
based on our lens models and evaluated at the cluster core 
radius. We see a good correlation and derive a best-fit 
slope of —0.6 ± 0.2. This compares well with the theoret- 
ical expected value from a tidal stripping model (Merritt 
1983) of -0.5: 
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Dark halos of the scale detected here indicate a high 
probability of galaxy encounters over a Hubble time within 
a rich cluster. However, since the internal velocity disper- 
sions of these cluster galaxies (< 250kms~^) are much 
smaller than their orbital velocities, these interactions are 
unlikely to lead to mergers, suggesting that the encounters 
of the kind simulated in the 'galaxy harassment' picture 
(Moore et al. 1996) are frequent and likely. In fact, high 
resolution cosmological N-body simulations of cluster for- 
mation and evolution (Ghigna et al. 1998; Moore et al. 
1996), find that the dominant interactions are between 
the global cluster tidal field and individual galaxies after 
z = 2. The cluster tidal field significantly tidally strips 
galaxy halos in the inner 0.5 Mpc and the radial extent 
of the surviving halos is a strong function of their dis- 
tance from the cluster center. Much of this modification 
is found to occur between z = 0.5-0. Detailed comparison 
of these results with tidal stripping of dark matter halos 
in cosmological N-body simulations will be presented in a 
forthcoming paper. 
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The prospects for extending this technique to larger 
scales within clusters in order to study the efficiency of halo 
stripping as a function of radius (variation of as a func- 
tion of radius) and morphological type are very promising 
with new instruments such as the Advanced Camera for 
Survey on HST. Multi-band imaging will enable photo- 
metric redshift determination for the background sources 



which will reduce one of the significant sources of noise for 
future analyses. 
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Fig. 1. — The results of the maximum-likelihood analysis for the HST cluster-lenses in our sample, <t* and correspond respectively to 
the central velocity dispersion and outer scale radius (identified as the tidally truncated radius) for a fiducial L* cluster galaxy in each of 
these clusters. The contours start at 1-cr and increase in l-cr increments from inside out in all 5 panels. The thick solid lines represent the 
contours of constant aperture mass and the best-fit values at marked by the cross. 




Fig. 2. — Scaling out the variation in the <t*'s (as per eqn. 4 above), the central density of the cluster evaluated at the core radius is 
plotted against the tidal radius. The errors plotted in rt are the 3-cr values. Performing a least-squares fit to the data points, the value of the 
best-fit power- law index t], where oc p^, is estimated to be —0.6 ± 0.2, (dashed line) in excellent agreement with theoretical expectations 
of ?7 = —0.5 (solid line). 
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